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Abstract. We present a general radiative transfer model which allows the Zeeman diagnostics of complex and unresolved 
solar magnetic fields. Present modeling techniques still rely to a large extent on a-priori assumptions about the geometry of 
the underlying magnetic field. In an effort to obtain a more flexible and unbiased approach we pursue a rigorous statistical 
description of the underlying atmosphere. Based on a Markov random field model the atmospheric structures are characterized 
in terms of probability densities and spatial correlations. This approach allows us to derive a stochastic transport equation for 
polarized light valid in a regime with an arbitrary fluctuating magnetic field on finite scales. One of the key ingredients of the 
derived stochastic transfer equation is the correlation length which provides an additional degree of freedom to the transport 
equation and can be used as a diagnostic parameter to estimate the characteristic length scale of the underlying magnetic field. 
It is shown that the stochastic transfer equation represents a natural extension of the (polarized) line formation under the micro- 
and macroturbulent assumption and contains both approaches as limiting cases. In particular, we show how in an inhomoge- 
neous atmosphere asymmetric Stokes profiles develop and that the correlation length directly controls the degree of asymmetry 
and net circular polarization (NCP). In a number of simple numerical model calculations we demonstrate the importance of a 
finite correlation length for the polarized line formation and its impact on the resulting Stokes line profiles. 
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1. Introduction 

The entire solar photosphere exhibits a rich structure of large 
and small scale magnetic features like sunspots, pores, facu- 
lae or plages. But except for sunspots and pores these mag- 
netic fields are not spatially resolved with present telescopes, 
although these fields clearly manifest themselves in high res- 
olution spectropolarimetric observations. With the improve- 
ment of spectropolarimetric sensitivity and spatial resolution 
over the last years it became clear that these unresolved mag- 
netic fields are much more ubiquitous than previously thought. 
This raises the question of the significance of these elusive 
and complex magnetic fields for the solar magnetism in gen- 
eral dSchriiver & Titldl2Q03b [Sanchez Almeidall2004l) and how 
these magnetic fields can be appropriately investigated by 
means of spectropolarimetric observations and Zeeman diag- 
nostics. 

The concept of magnetic flux tubes as the building blocks 
of solar surface magnetism has surely played a key role in 
our understanding of the unr esolved magnetic field of the so- 
lar photosphere (Stenflol 1976b . A wealth of diagnostic tech- 
niques based on the flux tube concept have been developed 
over the years and greatly improved our physical insight into 
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the nature of the photospheric magnetism (e.g. Solankil ll993l) . 
The interpretation of the Zeeman induced spectral line polar- 
ization in the context of the flux tube model rests on the idea 
of tube-like magnetic structures surrounded by an effectively 
non-magnetic atmosphere. A so called 1.5-dimensional radia- 
tive transfer model is then applied where a number of rays 
(line-of-sights) piercing through the 2- or 3-dimensional ge- 
ometry of the model atmos phere in order to obtain its spec- 
tral signature ( Stenflol 1994 ). From the standpoint of radiative 
transfer this approach provides a rather macroscopic treatment 
of the problem, as each line-of-sight possesses its own indi- 
vidual atmospheric structure and the averaging process for all 
line-of-sights (LOS) is performed after the actual integration of 
the transfer equation. Another problem here is, a relatively de- 
tailed knowledge about the number density and the geometrical 
arrangement of the flux tubes is needed in advance. 

A quite different approach to characterize the small scale 
and unresolved nature of the photospheric ma gnetic field was 
put forward by Sanchez Almeida et al. ( 19961) . In an attempt 
to reconcile the ubiquity of asymmetric Stokes V profiles with 
the underlying magnetic field structure, they postulated the mi- 
crostructured magnetic atmosphere (MISMA). In that model 
the degree of fluctuation and intermittency is much higher 
than in the classical flux tube picture and the magnetic field 
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is assumed to be structured over scales much smaller than the 
mean free path of photons. This approach could successfully 
reproduce many of the observed St okes V profile asymme- 
tries in quiet and active solar re gions dSanchez Alme ida 1997; 
ISanchez Almeida & LiteslfeoOOl) . In terms of radiative transfer 
modeling the MISMA approximation is a statistical approach 
which requires that the photons undergo a rapid and random 
fluctuation on smallest scales along their trajectory (line-of- 
sight). If the fluctuation rate per length scale is much higher 
than the (true) absorption rate this scenario represents the clas- 
sical microturbulent approach and allows a simplified treat- 
ment of the radiative transfer. The statistical averaging process 
(convolution) over all atmospheric structures can then be per- 
formed before the actual integration of the transfer equation. 
Despite its appealing simplicity in the way this approach treats 
the radiative transfer the idealized assumptions about the un- 
derlying atmosphere strongly limits the application of this ap- 
proach. As soon as only one of the ensemble structures exceeds 
the micro-scale criteria (and we will see later that this happens 
at surprisingly small scales) the microturbulent approximation 
fails. Moreover, magnetic microstructures can not in general 
be describe in the microturbulent limit. As elongated thin mag- 
netic structures close to equipartition have a preferential verti- 
cal alignment due to buoyancy forces the path length of an indi- 
vidual line-of-sight that traverse through such a magnetic struc- 
ture depends on the angle between the line-of-sight and the axis 
of the structure. The same structures that manifest themselves 
as microstructures in observations near the solar limb can be- 
come macrostructures in disk-center observations. This simply 
reflects the fact that the vertical and horizontal dimension of 
flux structures in the solar photosphere are very different. 

Over the last decade numerical magneto-convection simu- 
lations have reached a level of realism where they can provide 
valuable information ab out the three-dimensional structure of 
the photospheric plasma dStein & Nordlund|l998|;ISteiner et al 



1998; Schussl eHl200l IVogler & SchussldbOOa IVogler et al. 
2005; Schaff enberger et all 120051: IStein & Nordlundl l2006lh 



Recent magneto-convection simulations of plage like regions 
have shown that flux expulsion and convective field amplifica- 
tion may lead to structures in the supergranular n etwork which 



resem ble thin magnetic flux tubes or flux sheets ( Vog ler et al 



2005). For the quiet sun, numerical simulations indicate a lesser 
degree of orga nization with more fragmented and incoher- 
ent s tructures (iSchaffenberger et all 120051; IStein & Nordlund 
2006j). These simulations suggests that neither a predefined 
static macrostructured nor a pure microturbulent approach is 
an appropriate starting point for the interpretation of spec- 
tropolarimetric observations. The magnetic field seems to com- 
prise a broad range of structures which span from micro- (over 
meso-) to macro-scales. This is picture is supported by re- 
cent Zeeman and Hanle based diagnostics which are consis- 
tent with a continuous magnetic field strength and flux distribu- 



tion ( Khomenko et al . 2003; iDominguez Cerdena et al. 2003b; 
iTruiillo Bueno et al.ll2004 

So the present paper is an attempt to enhance the diagnostic 
capabilities of the Zeeman effect and to bridge the gap between 
the macro- and microstructured paradigm to obtain a more gen- 
eral mesostructured approach. Our approach is based on a sta- 



tistical description of the underlying atmosphere and its rele- 
vant parameters in terms of a Markov random field. One of the 
aims is to explicitly incorporate the spatial coherences and to 
account for the finite character of the underlying magnetic field. 
Although the model atmosphere is inherently non-plane paral- 
lel this approach allows us to formulate the transfer of polarized 
light with the help of a one dimensional stochastic transport 
equation. 

The idea of a stochastic radiative transfer approach is 
by no means new, several attempts were made to describe 
the line formation in r a ndom velocity fi e lds, for exam- 
ple lAuvergne et all dl973l): iGail et al] dl974l): IFrisch & Frisch 
dl976[) : iMagnanl dl985h : iGu et al.1 dl995): Nikoghossian etal. 
dl997l) The literature for polarized line formation in stochas- 
tic magnetic media is considerably shorter, there are only a 
few attempts known to the authors bvlDomke & Pavlovldl979l) : 
Landi DeglTnnocenti & Landolfi (2004) and very recently 
bv ICarroll & Staudd d2003l l2005albl l2006h and IFrisch etal 
( 20051 l200rj|) . The latter authors nicely demonstrated how the 
finite character of the underlying structures affect the line for- 
mation. But unlike most of the aforementioned authors we will 
not consider the solution under some limiting aspects, our in- 
tention is rather to derive a general statistical description of the 
magnetized atmosphere and then, using this description, to de- 
rive a (stochastic) transfer equation for polarized light. 

This paper is organized as follows: In Sect. [2] we intro- 
duce the general statistical concept of the mesostructured mag- 
netic atmosphere. A model is presented which explicitly takes 
into account the spatial coherency of the atmospheric param- 
eters by means of a Markov random field. Based on that de- 
scription, we derive in Sect.[3]the stochastic transport equation 
for polarized light. We present a formal solution in terms of 
a modified evolution operator and show that the microturbu- 
lent and the macroturbulent approximations are special limits 
of our more general mesostructured approach. In Sect. [3] we 
also show that the origin of asymmetric Stokes profiles and 
a net-circular-polarization (NCP) are natural consequences in 
a mesostructured magnetic atmosphere and that the degree of 
asymmetry depends directly on the characteristic length scale 
(correlation length) of the underlying magnetic structures. In 
Sect. |4] we give a brief description of the numerical realization 
of our stochastic approach. Then a number of numerical ex- 
periments follow, where we demonstrate how the arrangement 
of the underlying atmosphere and in particular the correlation 
length of the magnetic structures have a decisive impact on the 
line formation process and on the appearance of the resulting 
Stokes profiles. Sect.|5]finally gives a brief summary and draws 
the main conclusions. 

2. The Mesostructured Magnetic Atmosphere 

Due to solar convection we know that the photospheric plasma 
has a rather complicated structure in terms of its temperature 
and velocity distribution. This dynamic behavior also influ- 
ences the magnetic field, in particular in a regime where the 
plasma f3 is close to one. It is this complexity why we have 
chosen to pursue a rigorous statistical approach to describe the 
atmosphere and the underlying magnetic field. 



T.A. Carroll and M. Kopf: The meso-structured magnetic atmosphere 



3 



In the following we will present a discrete description of the 
atmosphere, following the concept of numerical magnetocon- 
vective simulations, we approximate the atmospheric structure 
by a discrete 3-dimensional lattice structure with a regular ar- 
rangement of nodes. To each node we will assign a number of 
physical parameters (a state vector) which describe the physical 
conditions at the particular locations of the nodes. Each node 
is moreover associated with a multivariate probability density 
function and a neighborhood relation. This neighborhood re- 
lation will be given in terms of a Markov model which then 
allows us to derive, for an arbitrary trajectory or line-of-sight a 
differential equation which describes the evolution of the prob- 
ability density function. This differential equation forms the 
basis for Sect. [3] where we describe the line formation as a 
stochastic process and derive the stochastic transport equation 
for polarized light. 

2.1. The atmosphere as a Markov random field 

Let us begin by introducing the state vector B for an arbitrary 
position (node) in the atmosphere as a vector quantity which 
includes all relevant atmospheric parameters, 

B* = {a\,a 2 ,...a n ) T . (1) 

The vector B shall contain atmospheric parameters like temper- 
ature, velocity, magnetic field etc.. We assume, that the param- 
eters a n are random variables, hence the vector B represents a 
multivariate random variable or a random vector respectively, 
for which we define a probability density function (pdf) of the 
form 

p(B) =p(a\, a 2 ,...a n ), (2) 

where the probability P(B < B < B + dB) of finding a realiza- 
tion of B within B < B < B + dB is 

P(B < B < B + dB) = p(a u . . .,a„) da\ . ..da„ 

= p(B) dB , (3) 

and which satisfy the following normalization condition: 

p(B) dB = 1 . (4) 

OO 

Here, we have assumed that B is a real-valued random vector 
and we want to keep this notation in the following without loss 
of generality. For the sake of brevity and clarity, we also will 
not distinguish between the stochastic variable and its realiza- 
tions (states) hereafter, furthermore the integration limits will 
be omitted as long as the integration is performed over the en- 
tire domain or state space respectively. 

As mentioned in the introduction, we will now consider a 
spatially discretized approximation of the atmosphere in terms 
of a three dimensional uniform grid structure. The grid, here- 
after called the random field, is made up by a number of grid 
points (nodes) where each node is associated with a random 
atmospheric vector B,, and its respective probability density 
function 

p Si {Bi) = p(Bj, s t ) = p(a h ,a h ,... a in , s,) . (5) 



The index s,- specifies the spatial position of the node. In or- 
der to obtain a general probabilistic description of the entire 
field, we need to take into account that the random vectors B, 
at their specific positions s, might be correlated among each 
other. The statistical state of the grid is therefore described in 
its most general form by a joint probability density for the en- 
tire field p(B\, Si\... ,B„, s„). This joint pdf or n-point distri- 
bution function can also be expressed in terms of the respective 
conditional probability density function 

p(Bi,si;...;B„,s„) = p(B\,S\ I B 2 , s 2 ; . . .; B„, s„) 

xp(B 2 ,s 2 ;...;B„,s„) , (6) 

which emphasizes the spatial dependency (correlation) of the 
field. The state vector B\ at s\ will in general depend on other 
state vectors within a certain region or radius around S\. In 
general it will be difficult to formulate an approximation for 
the joint probability density function or conditional probability 
density (O in terms of an appropriate n-point correlation func- 
tion. To reduce the complexity of the random field the simplest 
assumption is that of a completely uncorrected field. Thus, 
the realization of a state vector B, at the position is com- 
pletely independent and uncorrelated from all other realizations 
at other positions. With this assumption the joint probability 
density © factorizes into the products of the individual one- 
point probability densities of the respective state vectors, 

p(B u si; . . ;B n , s„) = p(B\,s\)p(B 2 , s 2 ) . ..p(B„, s„) 

= Y]p(Bi,sd- (7) 

i 

This description of the atmosphere in terms of a spatially un- 
correlated field where the atmospheric state vectors B, are ex- 
clusively described by their respective one-point pdf, would 
immediately lead us to a microturbulent approach. Our inten- 
tion here is, however, to go beyond this microturbulent concept 
in order to take - at least to the lowest order - correlation ef- 
fects into account. In the following we will therefore consider 
the particular type of a Markovian field structure. A Markov 
random field is characterized by a rather simple neighborhood 
relation which describes the spatial correlation of each node in 
the field with its neighboring nodes. The neighborhood N, at 
the position s, can be defined as follows 

Ni = [ Sj : || Si - sj \\< r) , (8) 

where r gives the spatial distance or radius between s, and the 
neighboring points sj. The Markov property of the random field 
can then be defined such that for each state vector B, at Si the 
following conditional probability density function holds 

p(Bi, Sj | Bj, sj) ; where i + j and sj e Ni . (9) 

The conditional probability of the state vector B, is thus only 
correlated to its neighboring states. The underlying assumption 
here is that all higher order spatial correlations are sufficiently 
small (and rapidly decay to zero) compared to first order cor- 
relation effects such that the atmospheric vector B, at Sj decou- 
ples from situations (states) which are further away than the 
immediate neighborhood N,. 
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2.2. A sample path through the Markov field 

In this work we will concentrate on the radiative transfer in 
a non-scattering atmosphere and therefore let us consider a 
particular trajectory - a sample path - through the above de- 
scribed random field. The sample path is described by a one- 
dimensional, unidirectional and contiguous sequence (path or 
trajectory) with st < s 2 < . ■ . < s„ through the random Markov 
field, where s, now determines the node position along our sam- 
ple path. Note, that the state space vector B can either be dis- 
crete or continuous, without loss of generality and for math- 
ematical convenience we assume in the following a continu- 
ous state space. Along the trajectory we encounter a sequence 
of realizations with B\, S\\ B2, s 2 \ B„, s„ where the random 
vector Bj in general adopts different values at different posi- 
tions St. This stochastic sequence along the trajectory defines 
a stochastic process where the state vector B has a distinct 
dependence on the spatial position s. From the Markov prop- 
erty of the random field we know that there exist a conditional 
probability density for the vector B, at each position s,- in the 
atmosphere which relates B, to the state vectors of the neigh- 
borhood. For the stochastic process along our sample path the 
three dimensional neighborhood at each position reduces to a 
one dimensional relation. In particular for our unidirectional 
path this means that the state B, at s,- is correlated only to its 
predecessor state B, \ at s ,„i . This represent s the typical prop- 
erty of a Markov process (Ivan Kampenlll992h and allows us to 
write for the conditional probability densities at s„ 

p(B n ,s n \B n -\,s n -\;...;B\,s{) = p(B„, s„ | B„_i, s„_i), (10) 

where s„ and s„_i are neighboring points along our sample 
path. The conditional probability which, in the framework 
of Markov processes, is also called transition probability ex- 
presses the property that the probability for a transition from 
the state B n -\ at s n -i to a state B„ at s„ only depends on the 
condition (state) at s n -\ and not on the (spatial) history prior to 
B„-\ . As can be shown, the followin g equation must hold for a 
Markov process ( Ivan Kampenlll992l) . 



p(B 3 ,s 3 \B uSl ) = I p(B 3 ,s 3 \B 2 ,s 2 ) 

xp(B 2 ,s 2 \B u si)dB 2 , (11) 

where the integration is performed over the entire state space. 
This is the Chapman-Kolmogorov equation which states that 
the transition probability between two states B\ and B 3 with 
S\ < S3 corresponds to the product of the transition probability 
between the initial state B\ and some intermediate state B 2 as 
well as the transition from this intermediate state to the final 
state B 3 , integrated over all possible intermediate states. 

In the following, we will concentrate on a homogeneous 
stochastic process. A homogeneous process is characterized by 
a transition probability p(B 2 , s 2 \ B\,S\) that only depends on 
the spatial distance As = s 2 - s\ between the two states and not 
on their particular spatial positions. We therefore omit the in- 
dices and consider, without loss of generality, s as a continuous 
parameter. 

Let us now define a transition probability which obeys the 
Chapman-Kolmogorov Eq. (fTTT ). To obtain such a model we 



have to make assumptions about the underlying (magnetic) 
structure of the atmosphere and how this structure determines 
the short scale behavior of the transition probability. Following 
the idea of thin magnetic flux tubes we can assume that the 
photosphere is characterized by sharp transitions between in- 
dividual magnetic structures and their non-magnetic surround- 
ing. This picture may be justified by the extreme small bound- 
ary layers of thin magne tic flux tubes wh ich are well below the 
photon mean free path ( Schussler Il986l) . The Markov process 
must therefore exhibit a distinct discontinuous character where 
abrupt jumps between different atmospheric regimes along a 
trajectory occur. A suitable definition of the transition probabil- 
ity for small As, between two states B' and B" can be obtained 
by the following discontinuous process 

p(B", s + As I B', s) = e - y(B ' )As 6{B' - B ") 

+ (1 -e-y (B ' )As )q(B"), (12) 

where y(B') is the fluctuation rate of the atmospheric struc- 
ture B' and q a prescibed spatial independent probability 
density function. This probability density, which has to sat- 
isfy the normalization condition ©, allows a prior character- 
ization of the field and transition properties. Note, that the 
prescribed probability density q is in general not the sta- 
tionary (spatially-invarian t) probability density of the process 
dBrissaud & Frisch 1974t) . The fluctuation rate can be ex- 
pressed with the help of the correlation length A to give 



y(B)=A-\B), 



(13) 



As can be seen from Eq. (fT2l the correlation length can be 
considered as the mean length scale of the structures. The pro- 
cess (fl2l), which is som etimes called Kubo-Anderson process 
( Frisc h & Frisch] 19761) describes how the probability of the 
state vector B" depends on the spatial distance As between B" 
and the initially known vector B' . The probability for staying 
in the initial state or atmospheric regime B' decreases expo- 
nentially with the distance As while the probability for a tran- 
sition to another atmospheric regime B" increases with As and 
is weighted by the given probability density q(B"). As the con- 
ditional probability density of the process depends on the spa- 
tial distance As we can use the Kubo-Anderson process ( TTZb 
and the Chapman-Kolmogorov Eq. (fTTT > to derive the following 
integro-differential equation (see AppendixlAl which describes 
the evolution of the probability density for B along a particular 
trajectory. 



ds 



dp(B, s) = j w(B ^ B , s) ^ ^ dB , 

- j w(B',B, s) p(B, s) dB' , (14) 



where w is the transition rate (transitions per unit length scale) 
and is defined by (IA.2I ). This equation which is also called 
the master equation of the process describes how the proba- 
bility density evolves along a given trajectory through the at- 
mosphere. The first term on the r.h.s. represents a source or 
gain term for the probability, which describes transitions from 
an arbitrary atmospheric state B' to the state B under consider- 
ation. The second term on the r.h.s. is a sink or loss term for the 
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probability density and describes transitions from the state B to 
other available atmospheric states. The master equation forms 
the basis for the following derivation of the stochastic transfer 
equation for polarized light. 

3. The stochastic transport equation for polarized 
light 

3.1. Derivation of the stochastic transport equation 

We now consider a particular photon trajectory along the line- 
of-sight in a non-scattering atmosphere, which will allow us to 
describe the radiative transfer in terms of a stochastic process. 
We begin by writing down the (deterministic) transport equa- 
tion for polarized light along a given ray path coordinate s in a 
plane parallel atmosphere for a frequency v, which is given by 



dl v _ _ 

— — - Ayly + Jy 

ds 



(15) 



where / is the 4-dimensional Stokes vector / = I,Q, U, V T 
whic h describes th e intensity and polarization state of the light 
(e.g. IStenfloll 1994b . For clarity, we will omit the frequency in- 
dices in the following. The absorption matrix K and the emis- 
sion vector j are given by 



K 



and 



( m iq nu nv 
t]q m p v -pu 

nu -nv m pq 
v m pu -pq m ) 



j = UiJoJuJvf = (m,nQ,nu,nv) T B ., 



(16) 



where T denotes the transpose of the emission vector and B 
stands for the Planck function. Here we have assumed that 
the emission originates under local thermodynamic equilibrium 
(LTE). For the standard case of a spectral line formed in the 
presence of a magnetic field, the detailed expression for the 
entries of the absorption matrix K a s well as fo r the emission 
vector j can be found for instance in iReesI (11987b . 

As our model atmosphere is described in statistical terms 
the observable quantity of interest in our case is now given by 
the mean value of the Stokes vector. We therefore define the 
mean Stokes vector by the first moment of the probability den- 
sity function p(I, B, s). Two things are noteworthy here, first, 
we have to include the spatial position s into the density func- 
tion since the value of the density function in general depends 
on the spatial position in the atmosphere since both / and B 
depend on s. Second, we have to include the Stokes vector it- 
self into the statistical description, because the Stokes vector 
is coupled to the atmospheric conditions via the transport Eq. 
< fl~5T > and therefore the Stokes vector itself becomes a stochastic 
variable which follows the same stochastic process as the at- 
mospheric vector B. The mean Stokes vector at the position a 
is then given by 



< I(s) > : 



II 



I p(I, B, a) dB dl. 



(17) 



We then separate the joint pdf p(I, B, s) by using the condi- 
tional pdf p(I | B, s) to write Eq. ( fT7] i as 



< I(s) > : 



/ p(I | B, a) p(B, s) dB dl . 



(18) 



The joint pdf can therefore be expressed as the product of the 
conditional pdf of the Stokes vector / which is conditioned 
on the atmospheric regime B at s and the structural pdf of 
the atmospheric vector B. We can evaluate the conditional pdf 
p(I | B, s) by realizing that this conditional pdf represents 
the deterministic solution of the transport Eq. ( fl3] l under the 
specific atmospheric conditions B. Thus expressing the condi- 
tional pdf with the help of the Dirac delta function we obtain 



p(I\B,s) = 6(1- KB, a)) = 6(1 -I(s)) 



(19) 



where Ks) is the formal solution of the initial value problem 
of Eq. $15[ given the atmospheric conditions B at s. The Dirac 
delta function of the Stokes vector has the following definition 

6(1 -I(s)) 

= 6(1 - I(s)) 6(Q - Q(s)) 6(U - U(s)) 6(V - V(s)) . (20) 
Eq. (fT8l can then be written as 

< I(s) >= j j 16(1- i(s)) p(B, s) dB dl , 



(21) 



Taking the derivative of ( fT8l with respect to the path variable s 
such that 



d<i> _ d 

ds ds 



s)) p(B, s) dB dl : 



(22) 



and interchanging the differentiation and integration operators 
we get 



8<I> 



-II 



I\—6(I-I(s))\p(B,s) 



+ /*(/- /(a)) — p(B, a) 

OS 



dB dl . (23) 



Taking further into account the differential form of the transport 
Eq. ( fl3b , we can write the derivative of the Dirac delta function 
d20T > with respect to the path variable s as 

^-6(1 - /(a)) = 2i£> V/ 6(1 - /(a)) 

OS OS 

= (-K B i + j B )V t 6(I-I(s)) 

= -(-K B i+j B )V I 6(I-i(s)) 

= -V/ [(-K B I + j B ) 6(1 - /(a))] , (24) 

where K B and j B are the absorption matrix and emission vector 
valid in the regime B, V/ is the nabla operator for the Stokes 
vector which is defined as 



1 

1 
nu 
jj_ 

ov 



(25) 
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Relation ( |24] |. allows us to write Eq. d23b as 

d<I> 

ds 

= /J -V'[(- K * I + hWI-foy)\p(.B,s)I dldB 



ff 



I S(I - I(s)) —p(B, s) dl dB. (26) 
ds 

If we now use the master Eq.dl~4b to express the spatial deriva- 
tive of p(B, s), Eq.(l26t can be written as, 

d <I > 

ds 

= J J-V/ [(-K B I + j B ) 6(1 -I(s))]p(B,s)I dldB 

+ ff{f w(B > B ^p( B '> s ) dB 'y 6(i ~i (s ^ didB 
J I [I w ( B '' B ^p( B > s)dB ') i6< ' i ~K s ^ didB - (2?) 

We then introduce the mean conditional Stokes vector (condi- 
tioned on the atmospheric regime B) which is defined as 



I p(I | B, s) dl . 



(28) 



The components of the mean conditional Stokes vector are 
given in analogy to the definition of the regular Stokes vector, 



I Y 1 \ 
1 B 
yQ 
1 B 

1 B 

yV 



(29) 



Taking into account Eq. ( fT~8b as well as the definition of the 
mean conditional Stokes vector we can integrate the first term 
on the r.h.s. of Eq. (|27| > by parts to obtain the following equation 

£ J" Y B p(B, s) dB = J (-K B Y B + j B )p(B, s) dB 

+ I [I W< ' B, B ) YB ' P(B '' S) dB ) dB 

w{B B)Yb P( B > s ) dB '^ dB > ( 3 °) 

where we have made the reasonable assumption that the atmo- 
spheric vector B is independent from the Stokes intensities. If 
we now neglect the outer integration over the state space B and 
divide Eq. ( f30b by p(B, s) as well as taking into account the 
definition of the transition rate (1A.21 i we obtain the following 
integro-differential equation for the mean conditional Stokes 
vector 

by r 

= -K B Y B +j B + j y(B')Y B 'p(B\ s) dB' 



I 



j(B)Y B q(B')dB'. (31) 



Here, we have introduced the modified fluctuation rate y(B) 
which is given by 

q(B) 



y(B') 



p(B, s)A(B') 



(32) 



Eq. ( f3Tb is the transport equation for the mean conditional 
Stokes vector Y B which can be solved by specifying the ini- 
tial conditions at the bottom of the atmosphere 



Y Bo (s ) = j 7 p (/o I B , so) dl 
= j I 6(I -I(s ))dlo, 



(33) 



as well as an initial probability distribution p(Bo, sq). Eq. OTt 
can be simplified under the assumption that the fluctuation rate 
y does not depend on the particular magnetic regime which 
then allows us to write 



dY B 

ds 



-KY B +j + y<I(s)>-yY E 



(34) 



The mean unconditional Stokes vector at the top of the atmo- 
sphere s, - our observable quantity - can then be readily ob- 
tained by a final integration of the mean conditional Stokes 
vector and the prescribed probability density q over the entire 
state space of B, 



< I(S,) > : 



(s t ) q(B) dB . 



(35) 



Eq. OTb describes the evolution of the mean conditional Stokes 
vector through the atmosphere from which the mean Stokes 
vector, the macroscopic (observable) quantity can be obtained. 
The stochastic transport equation for the mean conditional 
Stokes vector PD has the gener al form of a stochastic mas- 
ter equation (Ivan Kampen 1992) and describes the transport 
of photons which are subject to the familiar processes of ab- 
sorption (thermalization) and emission as well as two processes 
which describes the flow (fluctuation) between the different at- 
mospheric regimes or components of the atmosphere. To un- 
derstand the two new terms in the stochastic transport equation 
it is important to realize that Eq. OTb is conditioned on one par- 
ticular atmospheric regime B such that the integration process 
is entirely performed within that particular regime. Hence ,the 
two new terms in the stochastic transfer equation statistically 
couples the conditional transport equation for B to the other 
atmospheric regimes. In detail, the third term on the r.h.s. of 
Eq. (l3TT l gives the amount of intensity or photons which en- 
ters per unit length scale from all other atmospheric regimes 
into the regime B under consideration, this term therefore de- 
scribes an additional statistical inflow or source term, while 
the fourth term describes the loss of intensity per unit length 
scale due to transitions of photons to other possible regimes and 
therefore adds an additional statistical absorption. Due to the 
resemblance to the transport equation in a scattering medium 
both statistical terms can be considered as statistical scattering 
terms which describe the scattering of photons in and from the 
regime B. The stochastic transfer Eq. (l3TT l can be written in a 
more compact form if we combine the statistical source term 
with the emission vector j B and the statistical absorption with 
the absorption matrix K B . We can then define a modified (con- 
ditional) absorption matrix K B such as 



K B = K B + \y, 



(36) 
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where 1 is the unit 4x4 matrix. The modified (conditional) 
emission vector j can be defined as 



Jb = Jb + 7(B) Y B ' P {B',s)dB' . 



(37) 



This allows us to write the stochastic transport Eq. PTT ) in the 
following form 



9Yb - ky + 7 

-5 — - -H-B*B+jB 

OS 



(38) 



In the case of a discrete state space B we can obtain the dis- 
crete form of the stochastic transport Eq. PIT ) from the discrete 
master equation Eq. ( IA. 10b which gives 

^ = -KB,Y Bm +jB,„+J]y(B n )YB„p(B n ,s) 

n 

-Y j y(B m )Y Bn q(B H ). (39) 

n 

3.2. The formal solution 

A formal solution for the mean conditional Stokes 
vector can be obtained in analogy to the formal so- 
lution of the polarized radiative transfer equation 
dLandi DeglTnnocenti & Landi DeglTnnocenti Il985l) . Taking 
into account the definition of the modified absorption matrix 
(f36t the transfer equation (l38l in a purely absorbing medium 
can be written as 
dYn 

= -K B Y b • (40) 

os 

Since the statistical absorption is a symmetric operation and 
does not depend on the polarization state of the mean condi- 
tional Stokes vector we can introduce a modified conditional 
matrix attenuation or evolution operator which acts on the 
mean conditional Stokes vector Ob{s, s') at the position s' to 
give the mean conditional Stokes vector at the position s 

Y B (s) = d B (s,s')YB(s'), (41) 

and obeys the following limiting condition 

d B (s, s) = 1 . (42) 

Taking the derivative of d40b and taking into account the rela- 
tion ( |4TT i a differential equation for the conditional evolution 
operator can be written as 



dO B (s, s') 
ds 

Follo wing 



K B (s) O b (s, s') . 



(43) 



Landi DeglTnnocenti & Landi DeglTnnocenti 



(1985J), we can readily obtain the formal solution for the 
inhomogeneous case PIT ) by a direct substitution, such that we 
can write 



Y B (s) 



-f 

•J So 



Ob(s, s')j B (s')ds' + Ob(s, so)Yb(so) ■ 



(44) 



where /b(s') is the modified conditional emission vector ( ITTl i. 
In the special case of a constant modified absorption matrix Kb 
the conditional evolution operator can be written as 



O b (s, s') = exp \-(s - s')K B ] , 



where the exponential of the modified absorption matrix is 
given by its Taylor expansion. Thus we see that by intro- 
ducing a modified conditional evolution operator the for- 
mal solution for the stochastic polarized radiative trans- 
fer Eq. Pill is formally identical to the determin istic case 
(Landi De glTnnocenti & Landi DeglTnnocentill985 ). 



3.3. Micro- and macroturbulent limits 

In this section we show that the stochastic transfer equation 
OTb contains the microturbulent as well as the macroturbulent 
approach as limiting cases. 

We begin by considering the macroturbulent case which 
follows the idea that the turbulence or the stochastic character 
of the atmosphere is exclusively perpendicular to the line-of- 
sight. Any fluctuation along the line-of-sight is by far larger 
than the mean free path of the photons and is therefore negli- 
gible. On the other hand, the atmosphere is highly structured 
in the plane perpendicular to the line-of-sight. For observations 
with finite resolution, there always exists an ensemble of indi- 
vidual line-of-sights which traversing through different atmo- 
spheric regimes. The statistical variation within that finite res- 
olution element is usually described by a simple probability 
density function. Macroturbulence, whether for velocity field, 
magnetic fields or a more general atmospheric vector like B, 
certainly provides only a crude and rather simplistic represen- 
tation of a turbulent or fluctuating atmosphere with large coher- 
ent structures. In terms of our mesostructured approach, macro- 
turbulence can be described by a fluctuation rate that goes to 
zero or a correlation length that goes to infinity respectively. If 
we therefore consider Eq. QTl ) and take the limit as y — > or 
A — > 00 respectively, we obtain the transport equation for the 
mean conditional Stokes vector Yb in the macrostructured or 
macroturbulent limit 
8Yr 

lim — 2- = -KbYb+Jb- (46) 

,)->oo OS 

This leads to a deterministic transport equation where both sta- 
tistical scattering terms vanish and the conditional Stokes vec- 
tor in the regime B completely decouples from all other atmo- 
spheric components. A final integration over the state space at 
the top of the atmosphere (Eql35Tl gives the mean observable 
Stokes vector for each frequency 



< /(*) > 



Y B (s) q(B) dB . 



(47) 



In the case of a random macroscopic velocity field where B — v 
and a Gaussian LOS velocity distribution p(v), this equation 
turns into the familiar (macroturbulent) convolution integral. 

To derive the microturbulent limit, we have to take into ac- 
count that the transport equation ( IBTb was derived in the par- 
ticular limit of a small ratio of the path element A* to the cor- 
relation length A (see Appendix [A} and therefore the limit in 
equation ( l3TT l cannot be performed directly. Instead, let us first 
consider the underlying Kubo- Anderson process Eq.(fT2l with 
a depth independent fluctuation rate. If we take the limit as 
y — > 00 or A — ♦ we arrive at the relation 



(45) 



lim p(B", s + As \B',s) = q(B") . 

A — >0 



(48) 
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The transition probability from B' to B" is then exclusively 
determined by the probability of final state q(B") and, hence, 
entirely uncorrelated from the initial state B' . Note, that for a 
depth independent fluctuation rate the following relation holds 
Pstat - q, which states that the stationary solution of the mas- 
ter equation OTb is given by th e prescribed probability den- 
sity q dBrissaud & Frisch"lll974 ). The Kubo- Anderson process 
for this limit is therefore simply given by the overall station- 
ary (s-independent) probability density function p sla t{B) of the 
atmosphere. Since the Kubo- Anderson process is now indepen- 
dent of the spatial variable s it follows that the master equation 
(fT4l) is stationary as well and must be identical to zero. As the 
conditional Stokes vectors are in general Y B + this further 
implies that the sum of the second and third term on the r.h.s. 
of equation (TJTJ must be in an equilibrium state such that both 
will mutually cancel. This equilibrium can be expressed by the 
following relation, 



y 



J" Y B > PsmtiB') dB' = r J" Y B PstatiB') dB' . 



(49) 



From Eq. d48l ) and Eq. d32l ), we see that y — y and Eq. d49l can 
be written as 



y 



J" {Y B - - Y B 



) Pst JB')dB' = 0. 



From this relation it immediately follows that 
Y B =<1> . 



(50) 



(51) 



This relation is a direct consequence of the strong statistical 
coupling between the conditional Stokes vectors. If we insert 
the relation (IBTt into the transport equation OTb and taking into 
account relation d49l we can write the transport equation for the 
mean conditional Stokes vector in the microturbulent limit as 



dY B 

lim — — = -K B <I> +j B 

/!->°° OS 



(52) 



Note, that a strong statistical coupling of the conditional Stokes 
vectors, expressed by a small correlation length, leads to a rapid 
fluctuation of the atmospheric conditions along the photon tra- 
jectory. Due to this strong coupling and high transition proba- 
bilities, photons will fluctuate on very small scales many times 
from one atmospheric regime to another before they eventu- 
ally get thermalized or escape from the atmosphere. The high 
probability for a statistical absorption and scattering event com- 
pared to the thermalization probability provides an almost com- 
plete mixing on smallest scales such that the conditional Stokes 
vector will rapidly converge to the mean Stokes vector (see 
relation (BTTi). The microturbulent transport equation for the 
mean unconditional Stokes vector can then be derived from Eq. 
d35l l to give 



d<I > 

ds 



I 



K B <I> Pstat (B)dB+ j B p slat {B)dB. (53) 



I' 



using < I > = I, we can write the transport equation for the 
mean Stokes vector in the microturbulent limit finally as 



dl 

— — - < K > I + </>, 
ds J 



Note that, mathematically, the microturbulent limit is 
reached in the strict limit of A — > or y — > oo. However, as the 
complete mixing of the photons or the ratio of statistical scat- 
tering to true absorption is in fact controlled by the opacities in 
the individual atmospheric regimes the effective microturbulent 
limit will be reached for a finite correlation length (see Sect.|4}. 

3.4. Stokes profile asymmetries 

The asymmetry properties of Stokes profiles orig- 
inating in a stratified atmosphere were the sub- 
ject of several studies over the recent decades 



( Auer & Heasley 1978; Landi DegFInnocenti & Landolfi 
1983[ ISanchez Almeida et alJ Il989t ISolanki&Montavonl 



1 19931: iGrossmann-Doerth et all Food ISteineii I200C ; 
lLopez Aristd 2002 ). In the following, we show how the 
stochastic transfer equation for polarized light can be used 
to describe the generation of asymmetric Stokes profiles. 
We assume that the atmosphere comprises an ensemble of 
different magnetic regimes where each regime may have its 
own thermal properties and different LOS velocities. For 
simplicity, we assume that each regime or component is 
height-independent such that each conditional absorption 
matrix K B is constant with s. Note again, that the vector 
B c ontains all relevant atmospheric parameters . Accor ding 
to lLandi Pegl'Innocenti & Landi Pegl'Innocenti dl98lh for 
the Zeeman effect in an isolated triplet the seven entries 
of the absorption matrix possess clear symmetry properties 
with respect to the central line frequency or wavelength at 
any given depth. The point of symmetry of the absorption 
profiles in the absorption matrix are defined with respect to the 
central frequency Vo of the underlying Voigt- or Faraday-Voigt 
profiles. If we denote v as the distance from the the central 
frequency vo, then the entries of the absorption matrix K Vo 
satisfy the following relations 

7/M = m(-v) 

t]q(v) = t]q(-v) 

m(y) = nu(-v) 

rj v (v) = -j]v(-v) 
Pfi(v) = -Pq(-v) 
puiv) = -pu(-v) 

Pv(v) = Pv(-v) . (55) 

This is a direct consequence of the symmetry properties of 
the Voigt and Faraday-Voigt function. The absorption profiles 
r]j,r]Q,T] U ,p v are even functions regarding to the central fre- 
quency vo and the tjv,Pq,Pu, are odd functions. We recall 
that the evoluti on operator in the case of a constant absorp- 
tion matrix K dLandi DegFInnocenti & Landi Pegl'Innocenti! 
1985) can be written as 



0(s,s') = expt- K(s - s ' >] = J] 



-[K(s-s')r 



(56) 



n=0 



(54) 



Note that the evolution operator itself shares the same sym- 
metry properties about the central frequency vo as the entries 
of the absorption matrix K. We can now state, that any four 



T.A. Carroll and M. Kopf: The meso- structured magnetic atmosphere 



9 



component vector X whose components shares the following 
symmetry properties in the frequency domain 



X Y <* VIVO 

K x VVvo > 



(57) 



and which is subject to a transformation by means of the evolu- 
tion operator d56l > will preserve the symmetry properties. This 
can be verified by a direct multiplication, and is a direct con- 
sequence of the symmetry behavior (the even-odd relation) 
of the individual absorption profiles in the absorption matrix. 
Therefore, any evolution operator d56l > that acts on a vector 
X with the properties d57l) . and which shares the same central 
symmetry vo leaves the symmetry properties of the vector com- 
ponents unaltered. On the other hand, non-symmetric entries or 
entries with different points of symmetry in the evolution op- 
erator matrix will in general not preserve the symmetry of the 
vector X. 

For the sake of clarity we change the notation from a con- 
tinuous state space to a discrete, such that the individual atmo- 
spheric components B k can be better distinguished. If we now 
take the formal solution for the mean conditional Stokes vector 
Y Bk in a semi-infinite atmosphere we get 



Y Bt (s) 



f 



Bk (s, s')j Bk (s')ds' . 



(58) 



We then use the explicit form of the modified emission j Bk (s'), 
which is the sum of the thermal emission j T B herm and the inflow 
from the statistical scattering j B cat we can write Eq. d58l in the 
form 



jl h k erm (s)+f B c k at (s) 



j T B h k erm (s) + Y J 7(B k ,s)Y B ,p(B l ,s). 



(59) 



From the symmetry properties of the evolution operator 
Bk (s,s') and the thermal emission vector j J B " m (note, both 
acting in the same regime Bk) it is clear that a transformation of 
j J B * rm by the evolution operator Bk (s, •?') leaves the symme- 
try properties of j T B " m unaltered and therefore does not lead 
to any asymmetry of the conditional Stokes vector. This sim- 
ply reflects the fact that we have assumed that the absorption 
matrix for each atmospheric component is constant. We now 
concentrate on the statistical scattering term and replace Y Bl 
on the r.h.s. of Eq. ( 1591 by its formal solution and then insert 
the expression for the statistical scattering into Eq. d58l ) which 
gives 



%J — CO *J —CO i 



xj B ,(s")y(B k ,s')p(B,,s")ds"ds' . (60) 



Equation d60b leads to a sum of successive applications of the 
local evolution operators Bk and Bl . Despite the fact that 



each atmospheric component Bi has a constant absorption ma- 
trix, the individual atmospheric components differ in their at- 
mospheric parameters such that the matrix exponentials (evo- 
lution operators) in general do not commu te. We theref ore use 
the Baker-Campbell-Hausdorff formula dMagnusll 19541) to ex- 
press the product of two matrix exponential which involves all 
higher-order terms of the commutator, [A, B] = AB - BA, 

exp(A)exp(B) = exp(A + B + -[A,B] 

+^[A, [A,B]] - ±£B, [A, B]] + ...) . (61) 

Using this relation we can write Eq. d60b as 

Y Bk (s) = f f J] exp(K Bk (s' -s) + K Bl (s" - s') 

%J — CO %J — CO J 

+ l[K Bt ,K Bl ](s' -s)(s" -s') + ...) 

x }*,(*") y(B k , *') p(B,, s") ds"ds' . (62) 

Although each individual absorption matrix K Bl has its own 
central frequency v; (point of symmetry) due to the individ- 
ual atmospheric conditions, the superposition and multiplica- 
tion of the individual matrices will in general not retain or lead 
to a new symmetry in the compound matrix and therefore the 
resulting (combined) evolution operator will not possess any 
symmetry properties. This means that the modified emission 
vector j B ,{s") in Eq.(l62l. which again receives contributions 
from thermal emission and statistical scattering, is subject to a 
non-symmetric compound evolution operator of the form 

Btl (s, s') = ^ exp(K Bl (s' -s) + K Bt (s" - s') 



+ -[K Bk ,K Bl ](s'- S )( S "-s') + ...). (63) 

Regardless, of the symmetry properties of the statistical scat- 
tering vector j B cat (s") at s" - which itself is the result of the 
underlying mixed atmosphere - the symmetry of the thermal 
emission vector j T B herm (s"), which shares the local symmetry 
properties of the regime B\, is no longer preserved by a trans- 
formation with the compound evolution operator ( |63l and thus 
no general symmetry property for the conditional Stokes vec- 
tor Y Bk can evolve in this type of stochastic atmosphere. So 
we can conclude that it is the inherent asymmetry of the com- 
pound evolution operator (l63l that prevents any propagation of 
symmetric profiles. 

Yet another interesting fact that directly follows from Eq. 
(l62l is that the degree of asymmetry is proportional to the fluc- 
tuation rate y(Bt) or inversely proportional to the correlation 
length A(Bk)- This fact will be used in the following section to 
introduce the correlation length as a diagnostic parameter. 



4. Numerical experiments 

In this Section we perform a number of numerical experiments 
to demonstrate how mesostructures, structures of finite length, 
can affect the polarized line formation and the shape of Stokes 
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profiles. Particular emphasis in this section is given to the gen- 
eration of a asymmetric Stokes profiles and their dependence 
on the mean length scales of the underlying magnetic and non- 
magnetic structures. We show, how a wide range of structural 
length scales are not appropriately described by a microturbu- 
lent or macroturbulent approach and that significant deviations 
occur compared to our more realistic and accurate mesostruc- 
tured approach. 



The local absorption matrix /?, of the regime i is given by the 
following 4x4 matrix 





< 


-< 


< 






< 




-< 


Pi 




A 




< 


< 





Ri = 



And the local statistical scattering matrix 5, is given by 



(67) 



4.1. Numerical solution 

Before we begin with our model calculation we give a brief 
outline of the numerical procedure we have used to solve the 
stochastic transport equations. The numerical solution is based 
on the discrete form of transport equation for the mean condi- 
tional Stokes vector d39l which is given as 



dY Bi 



^ -K Bm Y Bm + j Bm + J] f(B n ) Y Bn p(B n , s) 

n 

-Y J 7{B m )Y Bm q{B H ). (64) 



Each of the transport equations for the mean conditional Stokes 
vector couples the 4 mean conditional Stokes parameter de- 
fined in Eq.d29li via the absorption matrix. Additionally, the 
transport equations are coupled via the statistical absorption 
and scattering terms to the n different atmospheric regimes. We 
therefore express the whole system by means of a An x An sys- 
tem matrix M which gives us the following system of differen- 
tial equations 



dY 

-r = MY + j , 

ds 



(65) 



where Y and / are the 4n-element Stokes vector and thermal 
emission vector. The system matrix M can be compactly writ- 
ten by introducing the local submatrices /?,- and 5,- to yield the 
following system of differential equations 



d 
ds 



YUs) 
TO 

TO 

YUs) 
TO 

to 
TO 



TO 

TO 
TO 
TO J 



R\ S2 S3 
Si R2 S3 

Si S2 /?3 



Si S2 S3 



s„ 
s„ 
s„ 



Rn 



yd 
yd 
yV 

Y l 

yd 
1 B, 
ylf 
1 B, 

1 B, 



Y' 

yd 

B„ 
ylJ 

B„ 
yV 

B„ ) 



1 /b, ^ 

J B\ 

L 



4 
i 



■(66) 



Si 















(68) 



The system is integ rated by a fourth-order Runge-Kutta method 
(|Press et al.l 1 1 9921) with appropriate boundary values for the 



mean conditional vectors Y Bm . The initial conditions are given 
according to Eq. J33I . 



Y Bm = S c (T iN )eo + 



dS c (r iN ) 



ds 



(69) 



where eo is the unity vector and S c is the Planck function in 
the LTE case. The lower boundaries are chosen for a sufficient 
large optical depth such that each of the different atmospheric 
regimes satisfy the condition r lA , > 2. 

To obtain the mean Stokes vector at the top of the atmo- 
sphere we finally sum up the individual contributions of the 
mean conditional Stokes vectors according to Eq. ( 1351 ). 

4.2. A mesostructured magnetic atmosphere - 2 
fluctuating components 

To demonstrate the importance of the correlation length for the 
line formation process we begin by considering a simple model 
of a stochastic atmosphere consisting of two different types of 
components. 

The arrangement of the different structures is such that we 
assume to have an ensemble of magnetic structures which are 
embedded in field-free regions. This picture is very similar to 
the often invoked two-component flux tube scenario often used 
for modeling unresolved magnetic structures in the solar net- 
work and internetwork. However, this scenario differs twofold, 
first, the resolution element comprises a whole ensemble of 
many individual structures and is not limited to a single static 
flux tube with a prescribed geometry, and second, the length 
scale of the individual magnetic and non-magnetic structures 
are finite. In the following we prefer therefore the term ensem- 
ble rather than components to emphasize the intrinsic complex- 
ity of the atmosphere. The line-of-sight length scale of the indi- 
vidual magnetic and non-magnetic structures in our model will 
be determined by the same correlation length L. 

The individual magnetic structures in our model share the 
same physical parameters, the same holds for the non-magnetic 
structures. Note also, that we make no particular assumptions 
about the geometry of the magnetic flux structures and the am- 
bient non-magnetic velocity field. For simplicity reasons we as- 
sume the same temperature stratification for the magnetic and 
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Fig. 1. The variation of the Stokes V asymmetry parameters (amplitude left; area right) as a function of the correlation length L. 
The dashed lines indicate the microturbulent limit for the asymmetry parameters calculated under the MISMA approximation. 
Note the rapid decrease of the Stokes V asymmetry with increasing correlation length. 



the field-free structures. This should be a reasonable choice for 
small scale magnetic structures (< 100 km in horizontal diam- 
eter) where thermal conduction by radiation efficiently smears 
out temperature differences. But even for larger structures it 
has no significant influence on the qualitative behavior of the 
here presented model calculations. The underlying tempera- 
ture and pressure struct ure is that of the quiet sun m odel atmo- 
sphere, HOLMUL, of iHolweger & Mueller! dl974l) . The mag- 
netic field strength within the magnetic structures is assumed 
to be 500 G (Gauss) with no internal LOS velocity. In the am- 
bient medium we assume a downflow LOS velocity of 550 m/s 
to mimic a intergranular downflow. The total filling factor of 
the magnetic structures, which corresponds to the probability 
value in our stochastic scenario, is set to 5%. Despite its sim- 
plicity this model should reflect the essential features of mag- 
netic flux structures embedded in intergranular lanes. For this 
numerical experiment we have used the iron line Fe I/16302A . 
The coupled discrete stochastic transport equation (f39b for each 
conditional Stokes vector is solved according to the numerical 
procedure described in Sect. ( 14. U . Note again, that we use the 
same correlation length L for both ensemble structures such 
that L is independent of the particular atmospheric regime. 

We now begin to change the correlation length L from very 
small (microturbulent) to very large values (macroturbulent) 
which results in a gradual increase in the (line-of-sight) order of 
the system (atmosphere). Such a decrease of the fluctuation rate 
should also affect the Stokes V profile asymmetry in character- 
istic way, because of the distinct dependence of the asymmetry 
on the correlation length L (see Sect. 13.41 ). To quantify the de- 
gree of Stokes V amplitude and area asymmetry we adopt the 
usual definition for the amplitude asymmetry 5a, 



5a — 



ml - \a r 



\a b \ + \a r \ 
and the area asymmetry 5A, 

L vdA 



5A = s b 



J A \V\dA 



(70) 



(71) 



where a b and a r are the extrema of the blue and red lobes of 
the Stokes V signal and s b the sign of the blue lobe. The inte- 
gration for the area asymmetry 5A is performed over the entire 
wavelength range A of the considered Stokes V signal. 

Figure [TJnow shows how the amplitude and area asymme- 
try of the Stokes V profile varies with the correlation length 
L of the magnetic and non-magnetic structures. For the broad 
range of meso-scales (20 - 1 000 km) we see the expected 
rapid decrease of the asymmetries with increasing correla- 
tion length L. This rapid decay already starts at a correla- 
tion length of approximately 15 km. For very small correla- 
tion lengths, the asymmetries (amplitude and area) saturate and 
begin to converge into the microturbulent limit. For a better 
comparison we have also calculated the expected amplitude 
and area asymmetries in the microtur bulent limit, where we 



have used the MISMA approximation (ISanchez Almeida et al 



1996) (marked by the dashed lines in Fig. [TJi. For very large 
correlation lengths (> 1 000 km), however, the asymmetries 
asymptotically converge to zero where they eventually reach 
- for a correlation length of approximately 10000 km - their 
macroscopic limits (complete symmetry). Fig. [TJ also reveals 
that unlike the mathematical limits, the microturbulent and 
macroturbulent limits are already reached numerically for fi- 
nite length scales. Figure [2] shows a sample of Stokes V pro- 
files from our model calculations and displays the variation of 
the Stokes V profile shape with the correlation length L. Note, 
in particular, the drastic change between the Stokes V profiles 
calculated for the L = 10 km and L = 100 km case. 



4.2.1 . About the microturbulent limit of photospheric 
structures 

Photospheric structures below 100 km are often considered to 
be optically thin and therefore have been treated in the micro- 
turbulent limit. From the latter section, however, we have seen 
that for a varying correlation length the shape, as well as the 
Stokes V profile asymmetries, undergo a rapid change, in par- 
ticular within the first 100 km. When there is a preferred inter- 
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Fig. 2. Sample of Stokes V profiles for a correlation length L 
= 10 km (upper left), L = 100 (upper right), L = 500 (lower 
left) and L = 1000 km (lower right). Note, the drastic changes 
in profile shape as well as in area and amplitude asymmetry. 

nal alignment e.g. by buoyancy forces the characteristic length 
scales of the magnetic structures in the horizontal and vertical 
can differ by more than an order of magnitude. This raises the 
question, to what extent structures below 100 km (measured 
along the LOS) can really be treated as micro turbulent in (po- 
larized) radiative transfer calculations ? 

From Fig.Q]we see that for a structural length scale of ap- 
proximately 15 km the atmosphere is effectively in a microtur- 
bulent state but apparent deviations from this microturbulent 
limit already occur from 20 km upward. For a characteristic 
length scale of just 100 km both asymmetries (area and ampli- 
tude) have lost already about 50 % of their initial (microtur- 
bulent) values. The rapid decrease of the area and amplitude 
asymmetries can be qualitatively understood from the inverse 
dependence of the asymmetries on the correlation length (see 
Sect. 13.41 1. The magnitude of the asymmetries is a direct con- 
sequence of the ratio of thermal absorption and statistical scat- 
tering. To quantify the degree of fluctuation and the statistical 
scattering occurring in the atmosphere in comparison to (reg- 
ular) absorption processes we define the statistical scattering 
probability as follows 



L seal — ! V 1 

K c + Kj + K s tat 

where the statistical scattering term K stat is the reciprocal of the 
correlation length L, k c the continuous absorption coefficient 
at line center and k/ the line absorption coefficient. We have 
calculated the statistical scattering probabilities for a number 
of different correlation lengths, where we have used the same 
atmospheric condition as in the previous section. In Fig. [3] we 
can see the run of the statistical scattering probability over the 
logarithm of the optical depth for different correlation lengths. 
For a better illustration and a better estimate of the line forming 
regions we have plotted in the background of Fig. |3]the wave- 
length integrated mean line depression contribution function 
for the circular polarization (solid line). The wavelength inte- 
grated mean line depression contribution function for the circu- 



Fig. 3. The statistical scattering probability for different corre- 
lation lengths over the optical depth. In the background, the 
wavelength integrated mean line depression contribution func- 
tion for the circular polarization is given (see explanation and 
definition in text). 

lar polarization is an integral measure that assess the contribut- 
ing layers of the emergent circular polarized signal in a stochas- 
tic atmosphere. The exact definition is given in AppendixlBl 

In very deep layers (log(Ts) > 0) of the atmosphere we see 
from Fig.[3]that for all correlation lengths but the smallest (10 
km, the effective microturbulent limit) the probability that pho- 
tons are statistically scattered into other regimes is very small. 
Most of the photons in these deep layers are subject to true 
absorption processes and will be thermalized before they can 
scatter into other atmospheric regimes. With increasing geo- 
metrical height or decreasing optical depth we see a gradual 
increase of the statistical scattering probability for all curves. 
But note that the increase of the scattering probability and the 
slope of the curve have a distinct dependence on the correla- 
tion length. In very high layers of the atmosphere the optical 
depth is so small that almost all the photons are subject only to 
statistical scattering regardless of their correlation length. Very 
interesting in that context are the statistical scattering proba- 
bilities of the L = 10 km (short dashed line) case and the L = 
100 km (dashed dotted line) case. The two slopes significantly 
differ from each other even though both may be regarded in 
classical terms as microturbulent. The L = 10 km curve for 
the effective microturbulent case already indicates a significant 
scattering probability for very deep layers, such that the entire 
line profile (wings and core) is affected by the strong statistical 
scattering. Whereas, the L = 100 km curve indicates a rather 
small scattering probability for the deepest line forming layers 
around Zog(Tg) = 0. At an optical depth where the integrated 
Stokes V contribution function reaches its maximum, the L = 
10 km curve has a value of approximately 70 % which means 
that photons are much more likely to encounter a transition to 
another atmospheric regime (70 %) than being thermalized in 
their current regime. Thus, most of the photons are statistically 
scattered many times before they are finally thermalized in a 
true absorption process. At the same optical depth the scatter- 
ing probability for the L = 100 km case is approximately 25 % 
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and therefore it has a less pronounced dependence on the sta- 
tistical scattering which results in a decrease of the line asym- 
metry of approximately 50 %. 

This behavior, even for a relatively moderate downfiow ve- 
locity of 550 m/s, as assumed here, demonstrates that the mi- 
croturbulent or microstructured (MISMA) approach does not 
give reliable results for atmospheric structures which exceeds 
the effective microturbulent limit. 

4.3. A magnetopause in a mesostructured atmosphere 

In the following, we want to take a closer look on the polarized 
line formation in an atmosphere that comprises a magnetopause 
like boundary. A typical example for such a scenario might 
be encountered in the vicinity of a rapidly expanding large 
flux tube structure ( Steinerl l2000). Similar structures with large 
overlaying canopies could be recently observed in magneto- 
hydrodynamic simulations dSchaffenberger et al]|2005l) where 
these canopies can cover large parts of the adjacent granular 
cell. The scenario we have assumed here for our model cal- 
culation is intended to be a crude approximation of a possible 
situation encountered in the photospheric internetwork. The ap- 
proximation we consider here can be conceived as an ensemble 
of flux tube like structures which are in pressure equilibrium 
with their field-free surroundings. These structures rapidly ex- 
pand until the individual structures finally merge together at a 
certain geometrical height to form a volume filling magnetic 
atmosphere. The salient feature of this model is the variation 
of the characteristic size of the individual magnetic structures 
with height such that these structures evolve from small to large 
scale structures. 

Again, we assume a simple 2-ensemble atmosphere with a 
number of magnetic structures which all shares the same phys- 
ical conditions as well as a number of non-magnetic structures. 
Furthermore, we assume that the magnetic structures are in 
pressure equilibrium with the non-magnetic surroundings, we 
therefore write the following relation for the height dependence 
of the magnetic field strength 



B( Z ) a B(z = 0) exp(-z/H B ) , 



(73) 



where Hb is the magnetic field strength scale height which is 
twice the gas pressure scale height and B(z = 0) the mag- 
netic field strength at the geometrical height z — 0. For Hb 
we have assumed a value of Hb = 250 km. The magnetic en- 
semble structures in our model scenario are mainly vertically 
orientated and we have adopted a magnetic field strength at the 
base z = of B(z = 0) = B(t = 1) = 500 G. For clarity 
and to focus on the effect of a varying correlation length, we 
have ignored any possible velocities in the magnetic and non- 
magnetic structures. Again, we have used the temperature strat- 
ification of the HOLMUL atmosphere for both ensemble struc- 
tures. The merging height of the canopies (the magnetopause) 
is assumed to be at a geometrical height of z ~ 450 km which 
corresponds to a logarithmic optical depth of log(r 5 ) w -3.65 
in the non-magnetic HOLMUL atmosphere. At this height the 
field strength in the magnetic structures has already dropped 
to a value of approximately 80 G. For the magnetic ensem- 
ble structures we assume a fixed probability value of 1 % at 



z — . The non-magnetic ensemble structures have a prob- 
ability value of 99 % at z = and we assume a fixed quasi 
macrostructured correlation length of L = 1000 km below the 
magnetopause for these structures. A large correlation length 
seems to be justified because of the much larger fill fraction of 
the non-magnetic structures. To mimic the rapid transition for 
the majority of the line-of-sights at the boundary of the mag- 
netopause from whereon the atmosphere is assumed to be fully 
magnetized, we set the probability value for the magnetic struc- 
tures to 100 % and their correlation length to L = 1000 km from 
log(Ts) » -3.65 on upward. 

For this scenario, we have calculated synthetic Stokes spec- 
tra for the two neutral iron lines Fe i /16301 A and Fe i /16302A . 
This is a prominent line pair which is often used for Zeeman 
diagnostics. Our particular interest here is to demonstrate how 
this prominent iron line pair depends on the length scale of the 
atmospheric structures below the magnetopause. For this rea- 
son we take a closer look on the differential behavior of the 
Stokes V profiles of the two iron lines. A good measure to as- 
sess the strength of the profiles is the amplitude ratio or the 
line-ratio method ( Stenflol 1994 ) which is commonly used to 
distinguish if the observed Stokes profiles are the result of an 
unresolved weak (sub-kG) or strong (kG) magnetic field. The 
amplitude ratio is defined as a 6301 / 6302 = a (V 6m )/a(V 6302 ) , 
where a(V) represents the amplitude of the respective Stokes 
V profile. Yet, another method going in the same direction to 
distinguish weak from strong magnetic field fields is the ratio 
of the longitudinal magnetic flux de nsity derived from the two 



iron lin es / 6301 / 6302 = Bf^/Bfj^ (iDominguez Cerdena et al 



2003al) . For the magnetic flu x density Bj nnl , we have us ed the 



following approximation (see lSocas-Navarro et al . 2004) 



B 



Long 



c JET y(/t) dA 

geff £°[dI(A)/dA]dA 



(74) 



where C is a constant with C = -1/(4.67 x 10 13 A 2 ) and g eff 
is the effective Lande factor of the line. The relation d74l ), the 



flux ratio f 



6301/6302 



as well as the amplitude ratio a 



6301/6302 



are 



direct consequences of the weak field approximation dStenflo 
1994 . 

Keeping all the atmospheric parameters fixed, except the 
correlation length of the magnetic structures below the magne- 
topause we have calculated the synthetic line profiles of the two 
iron lines for a series of correlation lengths. In Fig. [4] we have 
plotted the amplitude ratio a 6301 / 6302 and the longitudinal mag- 
netic flux density ratio y 6301 / 6302 ver the correlation length of 
the magnetic structures. At first glance we see a rather surpris- 
ing result, despite the fact that the magnetic field strength is 
intrinsically weak (B(z = 0) = 500 G) (and even rapidly de- 
creasing with height) the amplitude as well as the flux density 
ratio suggest for small correlation lengths a strong underlying 
magnetic field in the kilo-Gauss range. The ratio for the Stokes 
V amplitudes shows a value of 1 .0 for small correlation lengths, 
although we would expect for our intrinsically weak field that 
the amplitude ratio is similar to the ratio of the corresponding 
Lande factors o f the two iron lines which is approximately 0.67 
dSolankilll993l) . Moreover, one would not expect a dependence 
on the underlying length scale. As the magnetic flux density 
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Fig. 4. The amplitude ratio a 6301 / 6302 vs. the correlation length (left) and The ratio of the longitudinal magnetic flux densities 
correlation length (right). The amplitude ratio as well as the magnetic flux ratio exhibit a clear dependence on 
the underlying length scale of the magnetic structures. 



is also retrieved from the weak field approximation we also ex- 
pect for the flux density ratio a value close to unity regardless of 
the correlation length. Our results show that this is not the case. 
For a small correlation length the flux density ratio converges 
to a value of 1.41 while for a correlation length, L = 1000 km, 
a flux density ratio of 1.1 is obtained. It is obvious that with 
increasing correlation length the amplitude as well as the mag- 
netic flux density ratio decrease. For a correlation length of ap- 
proximately 300 km the value of the amplitude ratio and the 
flux density ratio converge to their (expected) weak field val- 
ues. On the other hand, we see that for structures which have 
correlation lengths smaller than 100 km significant deviations 
from the expected values occur. Once again, we see that the un- 
derlying correlation length - the characteristic length scale of 
the magnetic structures - has a decisive impact on the polarized 
radiative transfer and the resulting Stokes V profiles. To better 
illustrate the drastic profile variation we have plotted the Stokes 
V profiles for the two iron lines Fe i /16301 A and Fe i /16302A 
in Fig. [5] for a correlation length L = 10 km (upper row) and 
a correlation length L = 1000 km (lower row). This counterin- 
tuitive behavior of the amplitude and flux density ratios can be 
understood in terms of the MESMA approximation. Figure [6] 
shows the wavelength integrated mean line depression contri- 
bution function of the Stokes V profiles C« ; cp (see Appendix 
( IBl for the two iron lines. These stochastic contribution func- 
tions Cr-cp tell us where in the fluctuating atmosphere the main 
contributions to the Stokes V signals originate. We immedi- 
ately recognize the distinct shift in the contributing layers of the 
Stokes V signal. Figure[6]reveals that for a characteristic length 
scale of L = 10 km both iron lines receive their major contribu- 
tions predominantly from above the magnetic canopy which is 
located at log(Ts) m -3.65. Whereas, for a mean length scale of 
L = 1000 km the main contributions almost exclusively come 
from below the magnetopause. The reduced contribution in the 
lower layers for both iron lines in the L = 10 km scenario, 
which is even more pronounced for the Fe i /16302A line, is a 
direct consequence of the absorption of polarized line photons 
in the non-magnetic layers. The non-magnetic structures do not 





/V ] 









Fig. 5. The upper left and right display the Stokes V profiles 
for the Fe 1 16301 A (left) and Fe i ,16302 A (right) line calcu- 
lated with a correlation length of L = 10 km. Both profiles in 
the L = 10 km scenario have the same amplitude and therefore 
would falsely indicate the presence of strong magnetic fields. 
The lower profiles which show the case for L = 1000 km, ex- 
hibit the expected behavior of the profile amplitudes for an in- 
trinsically weak field. 



only fill a substantial horizontal fraction of 99 % of the resolu- 
tion volume, they also have, and this is even more important, a 
larger correlation length and therefore once photons are statisti- 
cally scattered into these field-free structures, they have only a 
very small chance of being backscattered into a magnetic struc- 
ture. Structures with large correlation length are only in weak 
statistical contact with their surrounding. Whether photons are 
being scattered into, or originate from a non-magnetic struc- 
ture, the likelihood for staying in the non-magnetic regime for 
the rest of their trajectory is much higher for the photons than 
making another transition to a magnetic structure. As the scat- 
tering probability is very high in the magnetic structures (for 
low correlation lengths), many of the photons carrying the po- 
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larized information are scattered into the non-magnetic regimes 
where a substantial number of them is absorbed due to the in- 
creased opacity in the non-magnetic structures. This absorp- 
tion of polarized photons in the non-magnetic structures is the 
main cause for the reduction of the contribution in lower layers. 
Beginning from the magnetopause the atmosphere is magnet- 
ically coherent (essentially macro-structured) with a probabil- 
ity value of 100 % and a correlation length of L = 1000 km. 
Although the contribution to the Stokes V signal from these 
high layers is very small in absolute terms, it is still significant 
compared to the strongly reduced contribution from below the 
magnetopause. Hence, the absorption below the magnetopause 
is so efficient that the small contributions from above the mag- 
netopause outweighs the contribution from below. 

The situation is quite different for magnetic structures with 
larger correlation lengths. Although the probability value of the 
magnetic structures below the magnetopause is still only 1% 
the enhancement of the correlation length leads to a decrease 
of the statistical scattering from the magnetic structures into 
the non-magnetic ones. Hence, less photons are being scattered 
into the denser non-magnetic structures and more polarized line 
photons are able to stay for their entire trajectory through the 
atmosphere in the same magnetic structure. The increasing co- 
herency of the individual magnetic structures leads to the more 
regular shaped form of the contribution functions (normal line 
formation) for the L = 1000 km case. But, what is the cause of 
the unusual line ratios for small correlation lengths ? Now, as 
we have understood why there is such a significant difference 
in the contributing layers with respect to the correlation length, 
let us turn to the question of what causes the small but deci- 
sive difference in the contribution functions for the two iron 
lines. From Fig. [6] we see that, although the two contribution 
functions for the two iron lines in the L = 10 km case show 
qualitatively the same behavior, there is a notable difference 
between both contribution functions. The contribution function 
for the Fe i /16301 A line (dashed line) shows a larger contribu- 
tion in the layer above the magnetopause than that of the Fe i 
/16302A line (solid line). As both lines in the L = 10 km case 
receive their main relative contributions from above the mag- 
netopause the formation height of the two iron lines play an 
even more important role than in the unstructured and coherent 
(e.g. macrostructured) case. The larger contribution of the Fe i 
/16301A line is exclusively due to its higher formation height 
compared to the Fe i /16302A line. This difference in the for- 
mation height which was also illustrated and exp lained in terms 
of response functions by iMartfnez et al.1 (120061) as well as the 



above explained pronounced importance of the contributions 
from the higher layers above the magnetopause, lead to the ef- 
fect that the Stokes V signal of the Fe i /16301 A line receives an 
absolute larger contribution than the Fe i /16302A line. As the 
characteristic length scale of the magnetic structures grows this 
effect loses its significance because the peak of the contribution 
function begins to shift downward below the magnetopause and 
reduces the overemphasized difference in the formation height. 
This demonstrates once more that the analysis of magnetic field 
strengths based on the iron line pair at 6300 A in a structured 
magnetic atmosphere must be regarded w i th cau tion, as was 
also pointed out recently by IMartfnez et al. I d2006l) . 



4.4. A quiet sun magnetic field strength distribution 

In the recent years the possible role of the so called in- 
ternetwork magnetic field has become increasingly appre- 
ciated. Despite its elusive character, the ubiquity of these 
weak flux fields could be observed in a number of high 
spatial and high sensi t ivity spectropolarimetr i c measurements 
dSigwarfh et all 1 19991: iLin & Rimmeld Il999t ISigwarthl |200l[ 
Lites|2002t Dom inguez Cerdena e t alJl2003atlKhomenko et 



x3 

tall 



2003; Socas- Navarro et alJ 12004). The significance of these 
fields for the solar atmosphere is not yet clear, but the 
mere fact that these fields cover most of the solar surface 
makes them probably an important ingredient for the struc- 
turing and dynamics of the higher layers of solar atmo- 
sphere dSchriiver & Title! 120031: ISanchez Almeida! |2004 . The 
characterization of the internetwork magnetic field in terms 
of empirical as well as statistical parameters from Hanle 
and Zee man effect measurements i s curre ntly under debate. 
Recently Dominguez Cerdena et al. (2006) proposed a set of 
empirical probability density functions which describe the dis- 
tribution of magnetic field strength in the internetwork. They 
emphasize that field strength in the kilo-Gauss range, even 
though they have low probabilities, play the dominant role in 
transporting most of the magnetic energy into the upper solar 
atmosphere. 

For empirical estimates of these pdf's from spectropolari- 
metric observations it is important to realize that line formation 
is an inherent three dimensional process and it is therefore cru- 
cial to take into account the horizontal as well as the vertical 
structuring of the atmosphere. In the following simple numeri- 
cal model calculation we want to us e the lognormal probability 



distribution function proposed by Do minguez Cerdena et al 



(2006) to demonstrate that the structural length scale of the in- 
dividual magnetic structures has an important impact on the 
resulting Stokes signal and that empirical probability density 
functions are strongly biased by the underlying model assump- 
tions. 

The following probability density for the magnetic field 
strength is used 



P(B) 



1 



yfn crB 



exp 



(In B- In B ) 2 



(75) 



The parameters Bo and <x are related to the first and sec- 
ond order moments respecti vely. We have adopted the value s 
cr - 1.7 and Bq — 38 G (see lDominguez Cerdena et al. I l2006h . 
The probability density function (l75T l for the magnetic field 
strength is shown in Fig. [7] note, that the probability for a field 
strength larger than 1000 G is only approximately 0.33 %. It 
is now crucial to realize that this empirical probability den- 
sity function comprises a broad range of different field strength 
regimes, a weak and probably purely turbulent range as sug- 
gested by Hanle effect measurements, a wide range of inter- 
mediate field strengths which are approximately in equipar- 
tition to their non-magnetic surroundings and a strong field 
range which is assumed to be populated by small intense kilo- 
Gauss structures. This picture is also supported by an increas- 
ing evidence for the existence of a small scale mixture of weak, 
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Fig. 6. The wavelength integrated mean line depression contribution function over the optical depth for the Fe i /16301 A (dashed 
lines) and Fe i /16302A line (solid lines) for a correlation length of L = 10 km (left) and for L = 1000 km (right). Note, the drastic 
changes of the contribution functions for both lines simply by varying the correlation length. For a correlation length of L = 10 
km the contributions to the Stokes V signal come predominantly from above the magnetopause. Note, the enhanced contribution 
for the Fe i /16301 A line in the L = 10km case (left). In the L = 1000 km case (right) one sees that for both iron lines the main 
contributions come from below the magnetopause. 




Fig. 7. Empirical probability density function for the magnetic 
field strength. 



intermediate and strong field s in a single resolution element 
dSocas-Navarro & Litesll2004 . 

But do these different field regimes all have the same char- 
acteristic length scale ? As magnetic buoyancy keeps stronger 
magnetic field structures more aligned to the vertical direc- 
tion than weaker field structures, the extent of the line-of-sight 
which traverses through an individual magnetic structure can- 
not be the same for all field strengths. There must be a de- 
pendence of the correlation length on the field strength. To 
demonstrate the importance of this effect, we assume the fol- 
lowing simplified scenario: for field structures below one kilo- 
Gauss field strength we adopt a correlation length of 10 km (ef- 
fectively microturbulent) and for structures with field strength 
larger than one kilo-Gauss we vary the correlation length in 
a range from 10 to 10000 km. Again, we have synthesized the 
two iron lines at /16300A and have used the temperature stratifi- 



cation of the HOLMUL atmosphere. To assess and quantify the 
influence of the correlation length on the Stokes V profiles, we 
have used the longitudinal magnetic flux density (l74l i. which 
gives us an estimate of the magnitude of the circular polarized 
signal. 

Figure [8] shows how the magnetic flux densities vary with 
the correlation length of the strong field structures. For both 
iron lines one can see from Fig. [8] the rapid increase of the 
flux density with increasing correlation length. Please note, 
that no other parameter but the correlation length of the kilo- 
Gauss structures has been changed in our scenario, the prob- 
ability density function in our calculation is for all correla- 
tion lengths the same. For very small structures (~ 10 km) the 
magnetic flux density saturates and converges into the micro- 
turbulent limit which corresponds to the value calculated for 
the probability density function d75l > under MISMA approxi- 
mation. For very large structures (> 1000 km) the value for 
the magnetic flux density converges into the macroturbulent 
limit. We see that just by increasing the correlation length of 
the strong field structures from micro-scales to macro-scales 
the magnetic flux densities (and, hence, the circular polariza- 
tions of both iron lines) are more than doubled. This, again, 
demonstrates how drastic the underlying length scale of the 
magnetic structures affects the process of polarized line for- 
mation. Although the underlying probability density function 
(|75T > has not changed, the effect of the strong field structures is 
more and more pronounced simply by increasing their correla- 
tion lengths. The growing correlation length is responsible for 
the additional magnetic flux derived from both spectral lines. 

Another intriguing aspect is the increasing difference of the 
flux densities calculated from both iron lines with increasing 
correlation length (Fig. [8}. The flux density for the iron line Fe 
i /16301A seems to grow faster than that of the Fe i /16302A 
line. This difference in the magnetic flux density ratio would 
suggest that the underlying probability density is dominated by 
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Fig. 8. The longitudinal magnetic flux density as a function of the strong field correlation length, derived from the iron line Fe i 
/16301A (left) and from the Fe i ^6302A line (right). 
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Fig. 9. The amplitude ratio between the two iron lines Fe i 
/16301A and Fe i /16302A as a function of the strong field 
correlation length. Note, that the underlying pdf is for all cor- 
relation lengths the same. 

the kilo-Gauss structures. This misleading effect is also high- 
lighted by Fig. [9] which shows the amplitude ratio r 6301 / 6302 
of the two iron lines as a function of the strong field corre- 
lation length. The amplitude ratio would also lead to the er- 
roneous conclusion that the underlying probability density is 
strongly biased towards kilo-Gauss structures, even though the 
kilo-Gauss structures make up less than 1% of the overall dis- 
tribution. 

So we see again the decisive role of the correlation length 
and that the simple derivation of empirical probability den- 
sity functions from spectropolarimetric measurements is not 
straight forward. There is no direct link to the number density 
of strong field components without accounting at the same time 
for the inherent length scale (correlation length) of the mag- 
netic structures. The magnitude and the profile form of the po- 
larized signal is therefore much more the result of the complex 
interplay between the horizontal probability density function 
and the characteristic length scale of the individual structures. 



5. Summary and Conclusions 

The ubiquity as well as the degree of asymmetry of observed 
Stokes V profiles are direct indicators of the inhomogeneous 
nature of the unresolved magnetic field in the solar photo- 
sphere. Stokes V profile asymmetries can be explained in 
terms of atmospheric gradients or more general in terms of 
the finite spatial extent (correlation length) of the underlying 
magnetic field. This finite character of the individual struc- 
tures is the essential starting point of the here presented con- 
cept of a mesostructured magnetic atmosphere (MESMA). In 
contrast to the proposed micr ostructured magnetic atmosphere 
(ISanchez Almeida et al . 1996) the MESMA approach does not 
rely on the assumption of a solely microturbulent magnetic 
field, the MESMA concept rather allows all kinds of different 
characteristic length scales to be present in the atmosphere. 

We have put forward a very generic statistical model (Sect. 
|2]) which is based on a Markov random field. The key idea here 
was to explicitly account for the spatial coherency of the un- 
derlying structures. This model facilitates the derivation of a 
stochastic transport equation for the conditional Stokes vector 
(Sect. [3]). By utilizing a modified evolution operator we could 
derive in analogy to the deterministic case a formal solution for 
the stochastic polarized transport equation (Sect. [3~2b . In Sect. 
I3.3l we have shown that the microstructured (MISMA) as well 
as the macrostructured approach are special limits of our more 
general mesostructured approach. Hence, the here outlined 
concept of a mesostructured magnetic atmosphere presents the 
natural extension of the MISMA approach. We could also show 
that asymmetric Stokes profiles are a natural consequence for 
an atmosphere where velocity and the magnetic fields are or- 
ganized on finite scales. Moreover, our stochastic approach al- 
lowed us to determine a direct relation between the character- 
istic length scale of the atmospheric structures and the degree 
of the asymmetry (Sect. [3~4l >. 

One of the most important parameter in the MESMA ap- 
proach is the correlation length (the mean length scale) of the 
underlying atmospheric structures. As the correlation length 
appears explicitly in the stochastic transfer equation for polar- 
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ized light it can be used as a diagnostic parameter. In a num- 
ber of numerical experiments we used exactly this parameter 
to demonstrate the importance of the correlation length for the 
polarized line formation. 

In a first numerical experiment (Sect. 14.2b we have ana- 
lyzed the line formation in a fluctuating atmosphere which is 
characterized by two kinds (ensembles) of different structures 
(magnetic and non-magnetic). Albeit this is a commonly in- 
voked scenario for the interpretation of unresolved magnetic 
structures our approach fully accounts for the lack of knowl- 
edge about the underlying organization and geometry. In these 
model calculations it could be demonstrated that the Stokes V 
profiles as well as the amplitude and area asymmetries have 
a clear functional dependence upon the correlation length of 
the underlying structures. We also showed (Sect. l4~2~TI ) that the 
effective microturbulent limit is already reached for a finite cor- 
relation length of approximately 15 km. On the other hand this 
means that for correlation lengths larger than 20 km the atmo- 
sphere can no longer be described in terms of a microstructured 
magnetic atmosphere. 

An atmosphere which is characterized by an extended mag- 
netopause like structure was investigated in Sect. d4.31 >. Again, 
our results demonstrate the distinct change of the Stokes V 
profiles and other derived quantities with varying correlation 
length. The line formation strongly depends on the structure 
below the magnetopause. The ratio of the Stokes V profile am- 
plitude and the ratio of the magnetic flux density for the two 
iron line Fe i ,16301 A and Fe i /I6302A are significantly altered 
just by changing the characteristic length scale of the magnetic 
field. These ratios falsely indicate for small correlation lengths 
the presence of strong kilo-Gauss magnetic fields, although the 
entire model atmosphere is intrinsically weak. Thus, the valid- 
ity and application of these ratios as indicators for strong mag- 
netic field structures is highly questionable in inhomogeneous 
atmospheres. 

In Sect. ( 14.4t we finally investigated the polarized spectral 
signature which originated from a field strength distribution 
given by a recently derived empirical probability density func- 
tions. We placed particular emphasis on the fact that the under- 
lying magnetic field cannot properly be described in the context 
of a field strength distribution without taking into account the 
characteristic length scale of the individual structures. As the 
field strength covers a broad range from very weak to strong 
fields, the characteristic length scale of these structures has to 
be different. If we realistically assume that strong field struc- 
tures possess intrinsic length scales that are larger than that of 
the weaker field structures we have shown that the strong field 
part of the distribution dominates in contributing to the result- 
ing Stokes V profile and would lead to an overestimation of 
the strong field structures. Therefore neglecting the correlation 
length (line-of-sight filling factor) of the underlying structures 
will lead to incorrect estimates for the magnetic field and flux 
distribution. 

The here introduced concept of a mesostructured magnetic 
atmosphere (MESMA) is a statistical attempt to cope with the 
underlying complexity of the magnetic field in the solar pho- 
tosphere in terms of polarized radiative transfer modeling. As 
the majority of the magnetic flux in the solar photosphere is 



still not resolved and the fundamental length scale of the mag- 
netic structures is unknown, the here presented stochastic ap- 
proach offers a viable and promising tool for the interpretation 
of spectropolarimetric observations and the diagnostics of pho- 
tospheric magnetic fields. 

First Stokes profile invers ions on the basis of the MESMA 
approximation were made by lCarroll & Sta ude (2005bi 120061) : 
ICarrolll(l2007h . They have demonstrated the feasibility of their 
approach which allowed them to estimate the characteristic 
length scale of internetwork and penumbral magnetic field 
structures. 
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Appendix A: The master equation for the 
mesostructured magnetic atmosphere 

In following we want to derive a differential equation (master 
equation) that governs the spatial evolution of the probability 
density along the spatial coordinate s. To derive such a master 
equation for the Kubo-Andreson process let us take the Taylor 
expansion of the transition probability ( fT2l at the coordinate s' 
which gives 

p(B", s' + As | B', s') = (1 - y(B')As) 6(B" - B ) + 

+ y(B')q(B")As' + 0(As' 2 ), (A.l) 

where we have assumed that y(B') is sufficiently small com- 
pared to As such that second and higher order terms of the ex- 
ponential are negligible. We then introduce the transition rate 
w(B",B') which describes the probability per unit length scale 
for a transition from an atmospheric regime B' to a regime B" . 
The transition rate can be obtained from the spatial derivative 
of the Kubo-Anderson Process (TLZb . hence, for a sufficiently 
small y(B) the transition rate directly follows from Eq. IA.ll 



w(B",B') = y(B')q(B") 



(A.2) 



To introduce the dynamics of the process, we use the Chapman- 
Kolmogorov equation (fTTT) . 

p(B", s' + As | B, s) 

p(B", s' + As | B', s') p{B', s' | B, s) dB' , (A.3) 



and replace the first factor under the integral on the r.h.s. with 
Eq. ( IA. 1 b to obtain 

p(B", s' + As | B, s) = (1 - y{B")As) p(B", s' \ B, s) 

+ As j w(B",B')p(B',s' | B,s)dB' , (A.4) 

Please note, that we have already performed the integration of 
the first term over the delta function with respect to B'. We 
rearrange Eq. dA.4l ) to write 

p(B", s' + As | B, s) - p(B", s' | B, s) 
As 

= -y(B")p(B" , s' | B, s) 

+ j w(B", B')p(B', s' | B, s) dB' . (A.5) 



The definition of the transition rate dA.2l ) allows us to write the 
fluctuation rate y{B") in the following form, 



y {B") = f w(B',B")dB' 



(A.6) 
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note here, that q in ( |A.2t satisfy the normalization condition the line depression contribution functio n is then formally iden- 
tic. This relation can now be used to write (IA.5I ) as tical to the non-stoc hastic case (see e.g. lSolanki & Brul si ll 994 



p(B", s' + As \B,s)- p(B", s' | B, s) 



Rees & Geers 1996). By using again the formal solution dB- It 



we can write for the line depression contribution function of 



As 

the mean conditional Stokes vector 
MB", B')p(B', s' | B, s) dB' _ _ ~ 

C BR (logr) = ln(10) t 0,(0, t) < 7 (Q) > ■ ^ 
w(B',B")p(B", s' | B, s) dB' . (A.7) < 

The wavelength integrated mean line depression contribution 
Taking the limit as As -> gives the differential equation function of the circular polarization can then finally be defined 
which describes the spatial evolution of the conditional proba- as 
bility density, 



dp(B", s' | B, s) 



logr) = r r 

J A JB 



\C BR ; V (B,A,T)\p(B)dBdA, (B.4) 



& s> where the integration of the inner integral is performed over 

= I w(B" B')p(B' s' | B s) dB' tne state s P ace & an d tne outer integration over the wavelength 

J domain A of the spectral line. 



- j w(B', B")p(B", s' | B, s) dB' . (A.8) 

To obtain the differential equation for the unconditional proba- 
bility density function we multiply ( IA.8b with p(B, s) followed 
by an integration over the entire state space B. After rearrang- 
ing the indices we finally yield the master equation for the 
Kubo-Anderson process, 



= j w(B,B')p(B',s') dB' 



dp(B, s) 
~~d~s 

- j w(B', B)p(B, s) dB' . (A.9) 

The discrete analogue of the master equation is given by 
dp(B, s) 



ds 



= YjW(B,B')p(B',s') 

B' 



-Y_ i w{B',B)p{B,s) , (A. 10) 



B' 

where the sum spans over the entire discrete state space of B' . 



Appendix B: The wavelength integrated mean line 
depression contribution function for circular 
polarized light 

The wavelength integrated mean line depression contribution 
function can be derived from the formal solution of the mean 
conditional Stokes vector d44l > which reads for a semi-infinite 
atmosphere 

lfl(0) = O B (0,T)j B (T)dT, (B.l) 

Jo 

where Jb(t) is the modified conditional line emission vector 
( ITTl i. From the formal solution we can define the contribution 
function to the emergent mean conditional Stokes vector as 

C B (logT) = ln(10)TO B (0,T)j B (r). (B.2) 

To derive the line depression contribution function of the con- 
ditional Stokes vector we introduce the line depression of the 
conditional Stokes vector as R = 1 - -j-^. The derivation of 



